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ABSTRACT 

We present observations of LI 155 and LI 148 in the Cepheus molecular cloud, taken 
using the Far Infrared Surveyor (FIS) instrument on the Akari satellite. We com- 
pare these data to submillimetre data taken using the Submillimetre Common-User 
Bolometer Array (SCUBA) camera on the James Clerk Maxwell Telescope (JCMT), 
and far-infrared data taken with the imaging photo-polarimeter (ISOPHOT) camera 
on board the Infrared Space Observatory (ISO) satellite. The Akari data cover a sim- 
ilar spectral window and are consistent with the ISO data. All of the data show a 
relation between the position of the peak of emission and the wavelength for the core 
of L1155. We interpret this as a temperature gradient. We fit modified blackbody 
curves to the spectral energy distributions at two positions in the core and see that 
the central core in L1155 (L1155C) is approximately 2 degrees warmer at one edge 
than it is in the centre. We consider a number of possible heating sources and conclude 
that the A6V star BD+67 1263 is the most likely candidate. This star is at a distance 
of 0.7 pc from the front of L1155C in the plane of the sky. We carry out radiative 
transfer modelling of the L1155C core including the effects from the nearby star. We 
find that we can generate a good fit to the observed data at all wavelengths, and 
demonstrate that the different morphologies of the core at different wavelengths can 
be explained by the observed 2 degree temperature gradient. The L1148 core exhibits 
a similar morphology to that of L1155C, and the data are also consistent with a tem- 
perature gradient across the core. In this case, the most likely heating source is the star 
BD197053. Our findings illustrate very clearly that the apparent observed morphology 
of a pre-stellar core can be highly dependent on the wavelength of the observation, 
and that temperature gradients must be taken into account before converting images 
into column density distributions. This is important to note when interpreting Akari 
and Spitzer data and will also be significant for Herschel data. 

Key words: stars: formation - stars: pre-main-sequence - ISM: clouds - ISM: dust, 
extinction - ISM individuakCepheus, L1155, L1148 



1 INTRODUCTION 

Stars are known to form in the molecular cloud cores 
that make up the densest parts of the interstellar medium 
(ISM). The initial conditions for the star-formation process 
are r e presented by pre-stellar cores (Ward- Tho mpson et al.l 
1 1994 I Andre. Ward-Thompson, fc Barsonvl |2000| ). A pre- 
stellar core is a gravitationally bound fragment of cloud that 
has not ye t formed a hydrostatic object (i.e. a protostar) at 
the centre (| Ward- Thompson et alJfeOQTt ). If we are to under- 
stand the processes by which these objects are formed, and 
the manner in which the go on to form stars, it is important 
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to obtain observational data showing the density structure 
of these initial conditions. 

Due to their location, embedded in the centres of 
molecular clouds, pre-stellar cores can only be observed in 
emission at far-infr ared (far-IR) and submillime tre (sub- 
mm) wavelengths (War d-Thompson et alj Il994n. and in 
absorption at mid-IR wavelengths ( Bacmann et aLlboOOj ). In 
the far-IR and sub-mm, we observe the thermal emission of 
cold dust grain s. Pre-stellar cores were first identified obser- 
vationally by IWard-Thompson et aL I <| 19941) in the sub-mm, 
using the James Clerk Maxwell telescope (JCMT). Sub-mm 
observations are ideal for observing pre-stellar cores, in that 
they are sensitive to the cold est, densest parts of molecula r 
clouds (for more details, see IWard-Thompson et al.ll2007n . 
In addition, these observations can be carried out from 
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the ground, which means that we can take advantage 
of the large collecting area and high angular resolution 
of large ground-based telescopes such as the JCMT. 
In this series of papers we are studying the detailed 
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However, we cannot accurately determine the tempera- 
ture of the emitting grains from sub-mm observations alone, 
as these observations all lie on the Ray leigh- Jeans slope of 
the spectral energy distribution (SED). This temperature 
information is crucial as it is necessary for converting the 
observed flux density into a mass distribution. In order to 
obtain this temperature information, we need observations 
shortward of the peak of the SED, i.e. in the far-infrared 
(A < 200/rni). These wavelengths are only obtainable from 
space-based telescopes. In this study we compare observa- 
tions of a pre-stellar core in the Cepheus Flare molecular 
cloud, using both ground-based sub-mm observations, and 
space-based far-IR observations from Akari. 

The Cepheus Flare covers a large region of the sky, and 
contains a number of local (< 500 pc) star-forming clouds, as 
well as more distant OB ass ociations. Ma ny of these clouds 
were originally classified by iLvndsl (|l962l ') based on the ob- 
scuration of background s tarlight. The Cepheus Flare was 
studied in more detail by lYonekura et alj (|l997l ) using the 
molecular tracer 13 CO to derive column densities, sizes and 
masses for the component clouds. Distance estimates for the 
nearby Cepheus Flare clouds range from ~ 200 to ~ 400 pc 
(iRacind Il968l; jViottil [l969l: ISnelll Il98ll: IStraizvs et al.lll992j ; 
Kun fc Prustilll993l; Ivan den Ancker et al.lll997l; [Kunlll998l ; 
Kun. Vinko. fc Szabadosll2000l ; iBalazs et al.ll2004r ), a spread 
which is comparable to the linear separation of the clouds 
on the plane of the sky. The Cepheus Flare region has been 
extensively reviewed bv lKun. Kiss, fc Ba log (2008). 

Together with the majority of star-forming clouds 
within 500 pc of the sun, the Cepheus Flare forms part of 
a larger structure called the Gould Belt, which is a ring of 
molecular clouds and OB associations, approximately 700 pc 
across. The Cepheus Flare clouds are located at Galactic 
longitudes 110 - 115°, and latitudes 10 - 20°. 

In this work, we study the L1155 and L1148 clouds 
in Cepheus using the Far Infrared Surveyor instrument on 
the Akari telescope. L1155 and L1148 are part of a sub- 
group called the L1147/L1158 complex, which also contains 
L1147, L1152, L1158, and the well studied outflow source, 
LI 157. The su b-group was determ ined to be at a distance 
of 325 ± 13pc l|Straizvs et al.lll992T ), based on a plot of ex- 
tinction vs. distance for a number of stars along the line of 
sight to the sub-group. 



2 OBSERVATIONS 
2.1 Akari Data 

Akari is a space-based telescope operatin g across the en- 
tire infrared band, from 2.4 /jm to 160 pm (|Murakami et al] 



l2007fl . The satellite was launched and is operated by the 
Japan Aerospace Exploration Agency (JAXA). Its primary 
mission is to make an all-sky survey with greater spatial res- 
olution and spec tral coverage than the da ta produced by the 
IRAS telescope dNeugebauer et al.ll 1984 ). The satellite is in 
a circular sub-synchronous polar orbit, allowing continuous 
scanning of the sky whilst avoiding pointing at the sun or 
the earth. In addition to the all-sky survey, a limited num- 
ber of 'observatory-type' pointed observations are possible, 
within these pointing constraints. 

The data described here were al l taken using the Fa r 
Infrared Surveyor (FIS) instrument l|Kawada et al.l [20071 ). 
This camera has four photometric bands: the N60 band is 
centred at 65 /im and has a bandwidth of 22 /im; the WIDE- 
S band is centred at 90 fim and has a bandwidth of 38 /im; 
the WIDE-L band is centred at 140 /im and has a bandwidth 
of 52 /im; and the N160 band is centred at 160 /im and has 
a bandwidth of 34 ^m. 

The N60, WIDE-S, WIDE-L and N160 detector arrays 
have 2x20, 3x20, 3x15 and 2x15 pixels respectively. The 
pixel scales are 26.8 arcsec for the two short- wavelength ar- 
rays, and 44.2 arcsec for the two long-wavelength arrays. 
The point spread function (PSF) at the four wavelengths 
are 37, 39, 58 and 61 arcsec respectively. The calibration ac- 
curacy is 20% for the N60 and WIDE-S bands, 30% for the 
WIDE-L band, and 40% for the N160 band. 

The FIS01 observing mode was used, which scans the 
arrays across the sky at an angle of 26.5° to the short axes 
of the arrays. The arrays are scanned at a relatively slow 
scan-speed of 15 arcsec s _1 (compared to the significantly 
faster all-sky survey observing mode). Each portion of the 
sky is mapped twice (known as a round-trip scan) in order 
to improve the detection redundancy. Each observation is 
made up of two round-trip scans, offset by half of the 480- 
arcsec field of view, in the cross-scan direction. The final 
map is made up of a mosaic of observations, each offset in 
the cross-scan direction by 450 arcsec in order to provide a 
uniformly sampled map. A reset interval of 1 s was selected 
to maximise the sensitivity, without saturating the detectors 
at the positions o f the b rightest sources. For more details, 
see iKawada et aL I (|2007D. 

The data were reduced using the AKARI official 
pipeline, version 20070714, developed by the AKARI data- 
reduction team (|Verdugo. Yamamura fc Pearsonll2007l ). 

2.2 SCUBA Data 

The sub-mm data presented in this study were obtained 
using the Submillimetre Com mon User Bolometer Array 
(SCUBA - iHolland et"ai]|l999l ) on the JCMT. SCUBA takes 
observations at 450 and 850 /im simultaneously through the 
use of a dichroic beam-splitter. The telescope has a resolu- 
tion of 8 arcsec at 450 fim and 14 arcsec at 850 /im. The data 
presented here were acquired from the JCMT data archive, 
operated by the Canadian Astronomy Data Centre. 

The observations were carried out over 5 separate nights 
between April 2002 and December 2003 using the scan-map 
observing mode. A scan-map is made by scanning the array 
across the sky, using a scan direction of 15.5° from the axis of 
the array in order to achieve Nyquist sampling. The array is 
rastered a cross the sky to build up a map several arcminutes 
in extent ijjenness fc Lightfootll200ol ). 
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Time-dependent variations in the sky emission were re- 
moved by chopping the secondary mirror at 7.8 Hz. The 
size of a scan-map is larger than the chop throw, therefore 
each source in the map appears as a positive and a negative 
feature. In order to remove this dual-beam function, each re- 
gion is mapped six times, us ing chop throw s of 30, 44 and 68 
arcsec in both RA and Dec |Emersonll 19951 ) . The dual-beam 
function is removed from each map in Fourier space by di- 
viding each map by the Fourier transform of the dual-beam 
function, which is a sinusoid. The multiple chop-throws al- 
low for cleaner removal of the dual beam function in Fourier 
space. The maps are then combined, weighting each map 
to minimise the noise introduced at the spatial frequencies 
that correspond to zeroes in the sinusoids. Finally the map 
is inverse Fourier transformed , at which point it no lon ger 
contains the negative sources jjenness fc Lightfootllioooh . 

The sub-mm zenith opacity at 450 and 850 /im was de- 
termined using the 'skydip' method and by comparison with 
polynomial fits to the 1.3 mm sky opacity data, measured 
at th e Caltech Submillimeter Observatory (|Archibald et al.l 
120021 ). The sky opacity at 850 /im varied from 0.18 to 0.37, 
with a median value of 0.26. These correspond to a 450 /im 
opacity range of 0.84 to 2.05, and a median value of 1.36. 

The data were reduced in the normal way us- 
ing the SCUBA User Reduction Facility, SURF 
l|jenness fc Lightfootl |2000| ). Noisy bolometers were re- 
moved by eye, and the baselines, caused by chopping 
onto sky with a different level of emission, were removed 
using the MEDIAN filter. Calibration was perf ormed using 
observations of the planets Uranus and Mars l|Orton et al.l 
1 19861 ; iGriffin fc Ortonlll993l ). and the secondary calibrator 
CRL618 !|Sandelllll994l ~ taken during each shift. We esti- 
mate that the calibration uncertainty is ±5% at 850 fj,m 
and ±15% at 450 fim, based on the consistency and 
reproducibility of the calibration from map to map. 



2.3 ISO Data 

The Infrared Space Observatory (ISO) satellite was a 
mid- and far-infrared telescope that was launched by 
the European S pace Agency (ESA ) in 1995 and oper- 



ated until 1998 (kessler et al.lll996l). The i maging photo- 



polarimeter (ISOPHOT- iKlaas et aLl Il994h was used in 
its over-sampled mapping mode (PHT32) to observe a 
number of pre-stellar cores, including L 1155C in Cepheus 
|Ward- Thompson. Andre, fc Kirkl I2002T ). Data were ob- 
tained at 90, 170 and 200 /J,m for this core with 
diffraction-limited resolutions o f 40 - 85 arcsec. See 
IWard -Thompson. Andre, fc Kirkl (|2002h for more details 
regarding the observations and reduction of these data. 
The absolute calibration errors were e stimated to be 30% 
l| Ward- Thompson. Andre, fc Kirkll2002l ). 

Th e characteristics of each of the fil ters used in 
Akari (IVerdugo. Yamamura fc Pearson! 120071), ISOPH OT 



3 RESULTS 
3.1 Akari 

Figures [TJ [2] and [3] show the L1155/L1148 region at 90, 140 
and 160 /im respectively. Wh ite contours showing the loca- 
tion of the optical extinction (jPobashi et al.ll2005h are over- 
laid on to each map. These contours are plotted to highlight 
any differences in the morphology of the Akari data at the 
three different wavelengths. The only source detected in the 
N60 data is the point source L1148-IRS. These data are 
therefore not shown here. 

There are a number of different sources see n in the data . 
These sources are labelled in Figure CD L1155 (|Lvndslll962l ) 
is located in the centre of the northern portion of the map, 
close to the northern extinction peak (Lynds clouds were 
discovered as opaque regions in the sky, therefore they would 
be expected to correlate with extinction peaks). LI 148 is 
located in the southern portion of the map, and is also seen 
as an extinction peak. 

Sin ce the study of extinction peaks in optical maps by 

iLvndsl (|l962t ). the L1155 clump has been subdivided us- 
ing C O and NH3 observations (e.g. iMvers. Linke. fc Benson! 
1983), and a number of these sub-clumps are seen in the 
Akari data. The extinction peak is centred on L1155C. 
L1155E is located 10 arcmin east of L1155C, and L1155D is 
approximately 12 arcmin southeast of L1155C. These have 
been labelled on Figure [1] 

We believe that L1155C is a pre-stellar core due to 
the detection of a centrally-condensed object at submil- 
limetre wavelengths with a high de nsity contrast with its 
surro undings (Section 13.21 - see also IWard- Thompson et aD 
2007). We also class this as isolated due to the large 
distance to the nearby cores, compared to more clustered 
star- forming regions su ch as Ophiuchus and Orion (e.g. 
i Motte. Andre, fc NerU Il998l; iNutter fc Ward- Thompson! 



ijKlaas et al 



given in Table [T] 



and SCUBA ijHolland et all Il999h are 



boOTl : ISimpson, Nutter, fc Ward-Thompsonl l2008). 

The other fragments of L1155 may also be pre-stellar in 
nature, or they may simply be transient starless cores. We 
have insufficient data to reliably determine their status, and 
so do not consider them in this paper. 

The emission at the far western edge of the 
map belongs to L1157, whic h contains the C l ass 
protostar IRAS 20386+675 1 (|Tafalla fc Bachillerl Il995l ; 
lLoonev. Tobin. fc Kwonll2007D . the driving source of a pow- 
erful bipolar outflow. The bright source seen on the north- 
eastern edge of t he maps i s the well-studied Herbig Ae/Be 
star PV Ceph (|Li et al.l 1 19941 ). which is the source of 
the parsec scale outflo w generating the HH 315 knots 
jArce fc Goodmanll2002l ). PV Ceph is no table for having an 
unus ually high velocity of over 20 kms -1 (|Goodman fc Area 
|2004| ). and is believed to have been ejected from the nearby 
star-forming cluster NGC 7023. 

L1148 contains the very low lum inosity Class (VeLLO) 
L1148-IRS (|Kauffmann et al.ll2005l ). This is a very low- mass 
Class protostar, and is believed to be either the precur- 
sor of a brown dwarf or a low-mass star, depending on the 
amount of gas it accretes in the future. This source is seen 
in both the N60 and the WIDE-S data. 

In Figures[T][2]and[3]we clearly see that the far-infrared 
emission peak of L1155C is not coincident with the extinc- 
tion peak. Instead, we see the far-infrared emission curv- 
ing around the south-west edge of the extinction peak. 
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R.A. (2000) 

Figure 1. Greyscale image of the 90 (im Akari WIDE-S data of L11 55 and L1148 with b lack contours overlaid. Contours start at 
9.7 MJysr -1 with intervals of 0.3 MJysr - 1 . White contours show the iDobashi et al. I <2005h visual extinction starting at Av=l with 
intervals of Av=0.5. The different clumps have been labeled. Note the offsets between the 90 fim peaks and the extinction peaks. The 
bright unresolved object near the extinction peak of L1148 is the very low luminosity Class (VeLLO) L1148-IRS. 




48 46 44 42 20:40 38 

R.A. (2000) 

Figure 2. Greyscale image of the 140 fim Akari WIDE-L data of L 1155 and L1148 with black contours overlaid. Contours start 
at 25 MJysr -1 with intervals of 5 MJysr -1 . White contours show the Dobashi et al. (2005) visual extinction starting at Av=l with 
intervals of Av=0.5. Note again the offsets between the 140 fim peaks and the extinction peaks. 



Akari mapping of pre- stellar cores 5 



1 o i=t- 




R.A. (2000) 

Figure 3. Greyscale image of the 160 fim Akari dat a of LI 155 and LI 148 with black contours overlaid. Contours start at 25 MJysr -1 
with intervals of 5 MJysr -1 . White contours show the lDobashi et all ll2005h visual extinction starting at Av=l with intervals of Av=0.5. 
Once again the far-ir emission is offset from the extinction peaks. 



The shortest wavelength shows this effect most prominently, 
peaking furthest away from the maximum extinction. The 
same is true in the L1148 cloud, where the far-infrared emis- 
sion wraps around the western edge of the extinction. The 
extinction peaks of the other sub-clumps of L1155 are unre- 
solved from the larger L1155C and L1148 peaks, therefore 
we cannot detect any offsets in these clouds. 

The far-infrared emission detected by Akari traces the 
density of the emitting dust weighted by its temperature, 
whereas the extinction map is based solely on the column 
density structure of the material providing the extinction. 
We therefore believe that the offset between the far-infrared 
and the extinction peaks can be explained as a temperature 
gradient across the core. This is considered in more detail 
in Section |4] 



3.2 SCUBA 

The 850 /im SCUBA data of L1155C are shown in FigureH 
Contours showing the 90 Akari WIDE-S data are in- 
cluded to highlight the position of the SCUBA core com- 
pared to the location of the core seen at other wavelengths. 
The sub-mm data shown have been smoothed from 14 arcsec 
to 60 arcsec. This was done in order to increase the signal 
to noise ratio of the data and to highlight the larger scale 
structures for comparison with the cloud models described 
in Section \4. 3 1 

We note that there is a small offse t between the peak o f 
the SCUBA data, and the peak of the lDobashi et all (|2005l 1 
extinction map. However, this offset is less than the 6 arcmin 



resolution of the extinction map, therefore we don't believe 
that this offset is significant. We also note that the morphol- 
ogy of the SCUBA data shows a very close agreement with a 
higher resolution extinction map generated from data taken 
with the IRAC camera on board the Spitzer space telescope 
(Chapman, 2009, priv. comm.). We therefore believe that 
the 850 /im peak represents the location of the maximum 
column density of L1155. 

None of the sources in the mapped area were detected 
at 450 ^m, therefore these data have not been shown here. 
The reason for the non-detection of any sources is proba- 
bly due to the poor atmospheric transmission at this wave- 
length. These data are used to provide upper-limits at this 
wavelength in Section U 

3.3 ISO 

Figure [5] sho ws the 90, 170 and 200 nm ISOPHOT d ata of 
L1155C from IWard- Thompson Andre, fc Kirkl (|2002l ), dis- 
played as contours overlaid on close-ups of the Akari 90 and 
160 fim maps. The 170 /jm ISOPHOT data are consistent 
with the 160 fim Akari data, showing a ridge of emission 
that stretches from the 160 /im Akari peak northwards. The 
200 fan ISOPHOT data peak at the northern extent of this 
ridge. The 90 /im ISOPHOT data miss most of the emission 
seen in the 90 /jm Akari data. The reason for this is two- 
fold. The 90 /im emission lies further south than the longer 
wavelength emission, and also the location of the detectors 
on the ISOPHOT focal plane meant that the 90 /im detector 
was offset, and in this case was unfortunately shifted slightly 
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Figure 5. The ISOPHOT maps of L1155C at (a) 90 fim, (b) 170 (im and (c) 200 /im, shown as thick black contours. The contour levels 
are (a) 18 MJySr -1 , with 0.5 MJySr -1 intervals, (b) 40 MJySr -1 , with 10 MJySr -1 intervals, and (c) 30 MJySr" 1 , with 10 MJySr -1 
intervals. The background greyscale and white contours show a close-up of the Akari maps at (a) 90 fim and (b) & (c) 160 (im, with 
contour levels as in Figures [T] and [3] The similarity between the ISO data and the Akari data is striking. 
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Figure 4. The SCUBA map of L1155C at 850 ftm with thin black 
contours overlaid. The data have been smoothed to a FWHM of 
60 arcsec to improve the signal to noise ratio and to highlight the 
larger scale structures in the data. The noise level in the central 
part of the smoothed map is 3 mJy/60-arcsec-beam. The thin 
contour levels are 3, 6 and 12-o\ We note that the far north and 
south portions of the map have lower integration times, therefore 
the noise-level is higher in these regions. The structure seen in 
these parts of the map are therefore less significant than in the 
map centre. The thick contours show the 90 fim Akari WIDE-S 
data, also smoothed to 60 arcsec for comparison, with contour 
levels as in Figure [T] 



to the northeast. However, the emission from the edge of the 
L1155C core detected by ISOPHOT at 90 fim is consistent 
with the Akari map at this wavelength. 

The ISOPHOT data therefore show the same trend for 
shorter wavelength emission to peak further away from the 
extinction peak, and back up the hypothesis of a tempera- 
ture gradient. The rms noise levels in each of the Akari, ISO 
and SCUBA maps is given in Table [1] 



Table 1. The characteristics of the different wavebands and the 
1-a rms in each of the maps. tlSOPHO T rms values taken from 
IWard-Thompson. Andre, fc Kirkl d2002l> 



Instrument 


Filter 


Central 


Bandwidth 


rms 




name 


wavelength 

w 


(H 


(MJy/ 


Akari 


N60 


65 


21.7 


0.5 




WIDE-S 


90 


37.9 


0.1 




WIDE-L 


140 


52.4 


0.7 




N160 


160 


34.1 


1.3 


ISOPHOT 


C-90 


90 


40 


0.7+ 




C-160 


170 


50 


5.7t 




C-200 


200 


30 


6. 5+ 


SCUBA 


450 


450 


66 


120 




850 


850 


22 


2.7 



4 ANALYSIS 

4.1 Temperature Gradient 

In the previous section we noted that the position of the 
emission peak of L1155C depends on the wavelength being 
observed, with the shorter wavelength radiation peaking fur- 
thest from the centre of the core. This trend is highlighted in 
Figure [6] which shows the position of the emission-peak as a 
function of wavelength. This shows that there is a monotonic 
shift of the peak with wavelength. 

As stated in the previous section, we believe that this is 
caused by a temperature gradient across the core. iKirk et al,l 
(2008) have shown that there are no embedded sources in 
L1155C, therefore the most likely cause of a temperature 
gradient is a heating source (presumably a star) located out- 
side L1155C. We discuss the possible identity of this heating 
source in Section \4. 21 

To investigate the temperature gradient hypothesis, we 
measured the flux densities at each of the wavelength^J at 

1 The 90 fim ISOPHOT data were not included in this analysis 
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Figure 7. SEDs for the centre and front positions. In both cases, the background flux density was removed from all points. The plotted 
line is a modified blackbody fit to the data. The error- bars show the calibration uncertainty at each wavelength (see Section [2$. 3-cr 
upper-limits are given for the Akari 65 (im and SCUBA 450 fim data. The best-fit temperatures are 11.9 K and 13.6 K for the centre 
and front positions respectively. See text for the parameters of each fit and for further details. 




1000 100 



Wavelength (jim) 

Figure 6. A graph showing the position of the peak of emission 
of L1155C as a function of wavelength, based on the Akari and 
SCUBA data. The position measured is at the closest point to a 
straight line joining the 90 and 850 (im peaks. The dashed line 
shows a fit to the data. The horizontal error-bars show the filter 
bandwidth. The vertical error-bars show the estimated positional 
accuracy in determining the peak of emission. 

two different positions. The first position is the centre of 
L1155C. To determine this position, we use the peak in the 
column density, which is the most likely the position of the 
volume density peak. We measure the position of the column 
density peak from the 850 /jm SCUBA map, because this 
wavelength is the least affected by the dust temperature, 
and is hence the most reliable tracer of column density. We 

because the data do not cover the positions of both the 850 /im 
SCUBA peak, and the 90 /mi Akari peak. See Section [3] for more 
details. 



note that there is material surrounding the peak, which is 
detected at 850 /im, and this material appears to contain 
sub-fragments. However, we do not include this material in 
our estimate of the position of the centre of L1155C, as we 
do not know what will happen to this material in the future. 
It could accrete onto the main L1155C core, fragment and 
form other stars, or it may simply disperse. 

As stated above, the morphology of the 850 /im emission 
shows excellent agreement with a higher resolution extinc- 
tion map, generated with data from the Spitzer space tele- 
scope. Therefore we don't consider t he offset between the 
850 /jm peak and the lower resolution iDobashi et alj (|2005h 
to be significant. The centre position is therefore taken to 
be 20:43:30, +67:52:51 (J2000). 

The second position (20:42:58, +67:46:27, J2000) is the 
peak of the 90 /im Akari map, which we hereafter refer to as 
the front of L1155C because if our hypothesis is correct, this 
region is the closest to the external heating source. The front 
position is located 0.5 pc away from the centre position. The 
flux densities at these positions were measured in circular 
apertures with a diameter of 2 arcminutes, and are given in 
Table [2 

The JCMT uses a chopping secondary mirror to remove 
the atmospheric emission (see Section [2.2[) . A consequence 
of this is that SCUBA does not measure the total power in 
the map, but only the power relative to the rest of the map. 
Akari and ISO, being satellites, have no need to chop and 
so they do measure the total power at each point on the 
sky. In order to compare the flux densities measured using 
ISO and Akari with those measured using SCUBA, all the 
fluxes were made relative to a common point, as far away 
from L1155C as possible. 

Spectral energy distributions (SEDs) were plotted for 
both the centre and the front positions. These SEDs are 
shown in Figure [7] Upper limits are shown for the 65 /im 
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Figure 8. A colour composite of the region around L1155 taken from the Digital Sky Survey (DSS). The red, green & blue channels 
show data taken using the infrared (850 nm), red (650 nm) and blue (480 nm) filters respectively. The white contours show the location 
of the Akari 90 fim emission (as for Figure [TJ- The positions of the candidate heating sources are labeled (see Table [3] for details). The 
front and centre positions are also marked with red circles. 



Akari data and the 450 fim SCUBA data, where no signal 
was detected. 

We fit a modified blackbody function to the data at 
each of the two positions in the norm al manner (see, e.g. 
I Ward- Thompson. Andre, fc Kirkf2002| ) . The best fits to the 
data are found to have a temperature of 11.9 K for the cen- 
tre position, and 13.6 K for the front position. The value 
of /3 (the dust emissivity index) used to produce these 
fits was 1.5. This differs from the canonical value of 2.0, 
which gave a poorer fit. However, we note that all values 
of f3 gave the same temperature difference of ~ 2 K when 
constrained to be the same at b oth p ositions. Following 
iKirk. Ward- Thompson, fc Andrei (|2007t) the critical wave - 
length (A c ) was set equal to 50 fj,m (see iHildcbrand 19831 ). 
Figure [7] indicates that the calibration may be marginally 
too low for the 160 /J,m N160 Akari data, as the flux density 
at this wavelength is systematically lower than the neigh- 
bouring wavelengths. 

These data are therefore in agreement with our hypoth- 



esis that the front of the L1155C core is warmer than the cen- 
tre. In the following section we consider the possible heating 
sources that could be the cause of this temperature gradient. 

4.2 Heating Sources 

The SIMBAD database was searched for all stars within 5 
degrees of L1155. This corresponds to 28 pc at the canoni- 
cal distance to the Cepheus molecular cloud of 325 pc. We 
only considered stars which were approximately the same 
distance away from the sun as the Cepheus molecular cloud 
(see Section [TJ, though we note that the uncertainty in the 
distances both to the cloud and to individual stars is large. 
Only stars of spectral type O, B, A and F were considered, 
as later-type stars would not be luminous enough to have a 
significant heating effect. The radiant flux due to each star 
at the position of the 90 /im intensity peak was calculated 
to determine the most likely candidates. Table [3] gives the 
details of these candidates. 
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Table 2. Background-subtracted flux densities measured at 
the centre (20:43:30, +67:52:51, J2000) and front (20:42:58, 
+67:46:27, J2000) positions in the Akari, ISOPHOT and SCUBA 
maps. Flux densities are measured in a 2 arcmin diameter aper- 
ture. 3-cr upper-limits are given for the Akari 65 fim and SCUBA 
450 fim data. 



Instrument 


Wavelength 


Flux density (Jy) 




(H 


Centre 


Front 


Akari 


65 


<0.4 


<0.4 




90 


0.4 


0.6 




140 


5.0 


5.3 




160 


3.4 


5.1 


ISOPHOT 


170 


8.0 


8.3 




200 


10.7 


9.5 


SCUBA 


450 


<96 


<96 




850 


2.7 


1.3 



The flux due to each star at the position of L1155C 
depends crucially on the distance between the star and the 
cloud. As this distance is not known accurately, upper and 
lower limits were estimated. The lower limit makes the as- 
sumption that the star and the cloud lie at the same distance 
from the sun, i.e. the sun-cloud-star makes an angle of 90°. 
It is assumed that this distance is the measured distance to 
the star, rather than the indirectly measured distance to the 
cloud. For the upper limit, we assumed that the sun-cloud- 
star makes an angle of 45° or 135°, therefore the distance 
between the star and the cloud is y/2 times larger than the 
distance in the plane of the sky. This assumption is based 
on the observed geometry of L1155C at the shorter wave- 
lengths. The core appears to being heated from one side, 
rather than the front or the back. 

Figure [8] s hows a colo ur composite of the digitized sky 
survey (DSS, Laskcr 199 4j), obtain ed using the SkyView in- 
terface (|McGlvnn fc Scollicklll994D . The red, green and blue 
channels represent the infrared, red and blue filters respec- 
tively. The contours of the 90 fim Akari data from Figure [T] 
are shown, to indicate the location of L1155C. The positions 
of the most likely heating sources are also highlighted. 

Table [3] indicates that the most likely heating source is 
the nearby A6V star BD+67 1263, as it produces an order 
of magnitude more heating flux that any other star. In the 
following section we use a radiative transfer simulation to 
determine if this or any other of the stars in Table [3] is 
appropriate for causing the temperature gradient that we 
see in the data. 



4.3 Radiative Transfer Modelling - L1155C 

The radiative transfer calculations were performed using 
the PHAETHON code, which is a 3D Monte Carlo radia- 
tive t ransfer code developed by [Stamatcllos & Whi tworthl 
l|2003h . PHAETHON simulates the heating of a cloud by in- 
jecting a large number of monochromatic luminosity packets 
(L-packets). These L-packets interact (by being absorbed, 
re-emitted or scattered) stochastically with the cloud. The 
injection point(s) of the L-packets represents the location of 
the radiation source(s) in the system. 

If an L-packet is absorbed, its energy raises the local 
temperature of the absorbing region. To ensure radiative 
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Figure 9. The azimuthally aver aged radial column-den sity profile 
of L1155C, determined from the Dobashi ct al. (2005) extinction 
data (open circles), with the best-fit Plummer-like sphere plotted 
as a solid line. See text for details. 



equilibrium, the L-packet is re-emitted immediately with a 
new frequency chosen from the difference between the local 
cell emissivity before and after the absorptio n of the packet 
jBiorkman fc Woodl l200ll: iBaes et all [2 005). For more de- 
tails, see IStamatellos fc Whitworthl l|2003l ). The model can 
produce images of the core at any wavelength, convolved to 
any resolution. 

The core is illuminated by a combination of the in- 
terstellar radiation field (ISRF), and the energy from the 
nearby heating source (see Table |3}. The ISRF is isotropic, 
whereas the energy from the nearby star is di rectional. For 
the ISRF we adopt a revised version of the IBlackl (|1994D 
interstellar radiation field (BISRF). The BISRF consists 
of radiation from giant stars and dwarfs, thermal emission 
from dust grains, cosmic background radiation, and mid- 
infrared emission f rom transiently heated small PAH grains 
I Andre et al.ll2003l ). The heating from the nearby star is pa- 
rameterised by the star's radius, surface temperature and its 
distance from the core. 

The initial shape of the core was taken to be a sym- 
metr i c sphere, with a plummer-like radial p rofile (|Plummerl 
1 19111 : IWhitworth fc Ward-ThompsonI 1200 ll ). This profile is 
approximately flat at the centre of the core and falls off with 
a power-law slope to large radii, according to the Plummer 
relation: 



n(H2)(r)=n(H 2 ) f iat 



flat 



( R2 fio.t + 



r 2\l/2 



(1) 



where r is the distance from the centre of the core, n{H2) fiat 
is the central density, Rfiat is the extent of the core in which 
the density is approximately uniform, and rj is the power- 
law slope at large values of r. A Plummer-like profile was se- 
lected because it provides a high-density inner region, which 
decreases to large radii, with a small number of parameters. 
This initial profile was used to model the density p r ofile o f 
pre-stellar cores by IWhitworth fc Ward-ThompsonI l|200ll) . 
We use the term Plummer-like because a true Plummer pro- 
file has rj — 5, whereas we treat this as a free parameter. 

The core is divided into a number of cells by spherical 
and conical surfaces. There are 70 spherical surfaces which 
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Table 3. Details of each of the possible heating sources in the vicinity of L1155. Column 2 shows the distance of each star from the sun. 
Columns 5 & 6 give the minimum and maximum distance between each star and L1155C (see text for det ails). Columns 7 fc 8 g ive the 
corresponding flux due to the star at the front of L1155C. Distance references: 1 . IStraizvs et al" | dl992r ); 2. iPerrvman et al.l dl997l ); 



Star 


R.A. 


Dec. 


Distance 


Spectral 


L 


Distance from front 


Flux at front of L1155C 




("20001 

V \J\J j 


(20001 


<oc) 


type 




of L1155C (pc) 


(10- 9 Wm- 
















Min 


Max 


Upper 


Lower 


BD+67 1263 


20:42:35.7 


67:36:35.1 


280 1 


A6V 


22 


0.7 


1.0 


1024 


521 


HD197053 


20:37:04.1 


67:30:15.5 


246 2 


B9V 


42 


2.7 


3.8 


148 


74 


SCKM 26 


20:40:28 


67:44:24 


480 1 


A5IV 


20 


1.8 


2.6 


148 


74 


BD+67 1265 


20:46:37.8 


67:58:58.4 


340 1 


A1V 


24 


2.4 


3.4 


108 


54 


HD198454 


20:46:32.2 


67:16:32.5 


344 2 


B9 


51 


.3.6 


5.1 


98 


49 


TYC 4447-2415-1 


20:40:31.2 


67:53:15.2 


200 1 


F6V 


2 


0.8 


1.2 


81 


41 


SCKM 48 


20:46:23 


67:57:24 


300 1 


Aim 


12 


1.9 


2.7 


78 


39 


HD198637 


20:47:40.2 


67:45:55.6 


330 1 


A3V 


19 


2.6 


3.7 


72 


36 


HD197794 


20:41:58.3 


66:54:43.7 


190 1 


A8IV 


12 


2.9 


4.1 


36 


18 


TYC 4447-614-1 


20:34:42.1 


67:56:54.3 


400 1 


A3IV 


31 


5.5 


7.7 


26 


13 


SCKM 2 


20:34:15 


67:47:54 


340 1 


A5V 


11 


4.8 


6.9 


11 


6 


HD194297 


20:20:42.3 


66:41:37.4 


360 2 




49 


15.0 


21.2 


5 


3 
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Figure 11. The best-fit model of L1155C 'observed' at (a) 90 fim and (b) 850 /im. The heating source is off the bottom of the figure. 
For comparison, the actual core is shown at (c) 90 fim and (d) 850 (im. These data have been rotated to match the observations and the 
model. The data have also been smoothed to a resolution of 60 arcsec to smooth out the fine details of the core that we do not attempt 
to model. In each case, the axes show the linear dimensions assuming a distance of 280 pc, with the origin at the centre position of the 
core. 
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Figure 10. The density profile of the best-fit model of L1155C. 
The scale on the right shows the logarithm of the volume number 
density, measured in units of cm -3 . The heating source is off the 
bottom of the figure. The centre and front positions are marked 
with a 'c' and an 'F respectively. 



are evenly spaced in radius, and 30 conical surfaces which 
are evenly spaced in polar angle. Hence the core is divided 
into 2100 cells. The number of cells used is chosen so that 
the density and temperature differences between adjacent 
cells are small. 

The value of Rfiat was determined from the 
iDobashi et all (|2005h extinction map of L1155C. Ideally, 
we would use the sub-mm data to estimate the extent of 
the core. However, the chopped nature of the sub-mm data 
means that they are only sensitive to spatial scales up to ap- 
proximately two-times the largest chop throw of 1 arcmin. 
This makes these data unsuitable for measuring the large- 
scale extent of the core. The extinction data however, are not 
based on chopped observations, and are also unaffected by 
the dust temperature. Therefore they trace the large-scale 
material well. 

An azimuthally averaged column density profile was 
first calculated from the extinction map. Equation [T] was 
then used to generate a volume density profile, which was 
decomposed to a column density profile under the assump- 
tion of spherical symmetry. Rfiat in equation [T] was varied 
to match the observed and derived column density profiles. 
The results of this are shown in Figure [9] for an Rfiat of 
0.35 pc. 

The central density (n(H2) fiat) was calculated from the 
peak flux-density of the SCUBA 850 /im data, which traces 
the densest region of the core, and is only minimally affected 
by the temperature g radient. The 6 arcmin resolution of the 
IDobashi et alj l|2005h extinction data make these data un- 
suitable for probing the very centre of the L1155C core, and 
hence for determining the central density. The uncertainty 
in the column density derived from the 850 /im data is ap- 
proximately ±75%, taking into account the uncertainties in 
the distance, dust temperature, and dust opacity. However, 
when this is converted to a volume density, the unknown 
physical size of the core along the line of sight introduces a 
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Figure 12. Graph of the emission from the best-fit model of 
L1155C at each observed wavelength, as a function of distance 
from the heating source, measured along a straight line connecting 
the core and the heating source. The flux densities at both the 
centre and the front positions are plotted. 



further uncertainty that is difficult to quantify. However, the 
aspect ratio of the core on the plane of the sky is approx- 
imately unity, so it is likely that the 3-dimensional aspect 
ratio is also of order unity. 

The parameters of the initial model, derived as de- 
scribed above, are Rfiat = 0.35 pc, and n{H2)fi a t = 2 x 10 4 
cm~ 3 . 77 was determined to be equal to 2, based on the fitting 
of the model to the data. 

However this model was unable to produce a good fit 
to the data. To improve the fit, the density profile was 
skewed towards the front of the core (i.e. towards the heating 
source). The form of this non-spherical Plummer-like profile 
is given by: 



A 



n(H 2 )(r,8) =n(H 2 )fiaf 



i*(6/2) 



1 + 



( Rfiat ) 



(2) 



(Stamatellos et al. 2004). The density profile for this model 
is shown in Figure [lOl The azimuthal angle 8 is defined with 
its origin pointing away from the heating source (the +y 
direction in Figure [10}. The parameter A determines the 
"south-to- north" optical depth ratio e, i.e. the maximum 
optical depth from the centre to the surface of the core 
(which occurs at = 180°), divided by the minimum op- 
tical depth from the centre to the surface of the core (which 
occurs at 8 = 0°). The parameter p determines how rapidly 
the optical depth rises with increasing 8, from the north at 
8 = 0° to the south at 8 = 180° . For the best-fit simulation 
we use Rfi a t = 0.35 pc, n(H 2 )fiat = 2 x 10 4 cm -3 , p = 4, 
A = 56.07, e = 1.5 and r\ = 2. Due to these additional free 
parameters, we accept that this model is not a unique solu- 
tion. However, the amount of skewing required is a relatively 
minor perturbation on an accepted basic pre -stellar density 
profile (|Whitworth fc Ward-Thompson|[200ll ). 

We also note, that this model is designed to reproduce 
the large-scale appearance of the L1155C core at each of 
the observed wavelengths. The model does not attempt to 
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Figure 13. The temperature of the best-fit model of L1155C, 
plotted as a function of the distance from the centre of the core 
in the direction of the heating source. Note the increase in tem- 
perature towards the front of the core. 



model any of the substructure of L1155C that is apparent 
in, for example, the 850 data. 

The results of the radiative transfer modelling are 
shown in Figures 1111 1121 and 1131 Figure [TJJ shows the 
model core 'observed' at (a) 90 and (b) 850 ^m using the 
PHAETHON code. For comparison, the data for the same 
two wavelengths are shown in Figure [TTJ (c) and (d) . These 
data have been rotated to highlight the similarities with the 
model, and a qualitative agreement between the model and 
the data at both wavelengths is clear. 

Figure [12] quantifies the similarity between the model 
and the observations by taking a 1-D cut through the model 
core at each wavelength. Each line shows the emission from 
the model core as a function of distance from the heating 
source, measured along a line connecting the core and the 
heating source. 

The flux densities at the centre and the front positions 
(listed in Table[2]| are also plotted on this figure for compari- 
son. There is very good consistency between the observations 
and the model at all wavelengths. 

The temperature profile of the best-fit model is shown 
in Figure [13] The temperature rises from ~ 10 K in the cen- 
tre of the core to ~13 K at the edge of the core nearest to 
the heating source. These temperatures are lower than we 
measured using the SED fits, especially at the central posi- 
tion of the core. However, this is expected in an externally 
heated core because the observations sample warmer mate- 
rial in front of the core as well as the colder material at th e 
centre |Stamatellos. Whitworth. fc Ward- Thompson|[20Q7h . 

We note that a change in the properties of the emitting 
dust grains might also affect the appearance of the core at 
different wavelengths, but the shift in grain properties would 
have to be dramatic, and a modest temperature variation 
seems far more likely. 

To conclude our study of the L1155C core, our hypothe- 
sis of a temperature gradient causing the observed change in 
core morphology at different wavelengths is confirmed by the 
SED fitting of the data at two positions in the core. In addi- 
tion, a radiative transfer model of L1155C with a geometry 



based on observed parameters shows that this temperature 
gradient can realistically be caused by a nearby bright star. 
This model successfully reproduces both the observed core 
morphology and measured flux density at each wavelength, 
and also the observed temperature gradient. 



4.4 Radiative Transfer Modelling - L1148 

The L1148 core contains the very low luminosity Class 
(VeLLO) L1148-IRS, and therefore represents an interme- 
diate phase between a pre-stellar core, and a protostellar 
core. The protostar has formed, but its luminosity is not yet 
strongly influencing the core. 

The L1148 core exhibits a similar morphology to the 
L1155C core in the Akari data, with the far-infrared emission 
wrapping around the peak of the extinction map (see Figures 
[TJ [2] and [3]) . We hypothesise that this is caused by the same 
process as in L1155C and that the observed morphology in 
the far-infrared is being affected by a temperature gradient. 

The data that we have available to constrain the model 
of this core are the 90, 140 and 160 (im Akari data over 
the whole core, and 850 /im data of the centre of the core 
only. The 850 fim data for the centre of LI 148 are shown in 
Figure I14H . and show that the core is centrally condensed, 
and peaks somewhat behind the crescent shaped 90 /im peak 
(see Figure PHf ). as is the case in L1155C. 

These data allow us to fit a modified blackbody curve 
to the SED at the centre position of L1148, as we did for 
L1155C in Section l"4.1l The best-fit blackbody has a temper- 
ature of 11.5 K, similar to the centre of L1155C. However, 
the lack of a long wavelength data-point for the front posi- 
tion means that we cannot constrain the blackbody at this 
position. Hence we cannot establish the presence or absence 
of a temperature gradient in LI 148 other than by analogy 
of the similar morphologies of L1148 and L1155C at far- 
infrared and sub-mm wavelengths. 

Figure [5] shows that the star BD+67 1263 is not in a 
suitable location to heat the L1148 core where we see it most 
brightly at 90 fim. A more likely heating source is the star 
HD197053 (see Table [3] for details). This star is located to 
the west of the L1148 core and has a high heating flux of 
between 300 and 600xl0~ 9 Wm 2 at the front of the core. 

We modelled the emission from the L1148 core with 
the PHAETHON radiative transfer code, using the same 
geometry as the best-fit model for L1155C. Again, the extent 
of the core and the p eak density were determined using the 
iDobashi et al.l (|2005l ) extinction data, and the 850 /im data 
respectively. 

The results of this model are shown in Figures [TJ] and 
1151 Figure [14] shows that there is a qualitative agreement 
between the model and the data. As for L1155C, the model 
reproduces the centrally condensed emission at 850 /jm, and 
the region of strong emission at the front of the core at 
90 /jm. The temperature of the model core increases from 
~ 10 K at the centre, to ~ 13 K at the front of the core. 
Again, the central temperature in the model is expected to 
be lower than the measured value. 

The agreement between the model and the data is 
shown quantitatively in Figure [TS] This shows the emission 
from the core as a function of distance from the centre in the 
direction of the heating source. The model is compared to 
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Figure 14. The model of L1148 'observed' at (a) 90 /im and (b) 850 fim. The heating source is off the bottom of the figure. For comparison, 
the actual core is shown at (c) 90 fim and (d) 850 fim. As for Figure [TT1 the data have been rotated to match the observations and the 
model, and smoothed to a resolution of 60 arcsec. In each case, the axes show the linear dimensions assuming a distance of 250 pc, with 
the origin at the centre position of the core. 



the measured flux densities at the centre and front positions 
of L1148, and shows a good agreement. 

In conclusion, though we have fewer data for the L1148 
core, the data are consistent with the core being externally 
heated. This heating is dramatically influencing the appear- 
ance of the core at far-infrared wavelengths. This is espe- 
cially interesting, because LI 148 harbours the ver y low lumi- 
nosity Class L1148-IRS (|Kauffmann et al.ll2005T ). However, 
the data are consistent with the core being heated externally 
rather than there being a temperature gradient caused by 
this heating source. This is possibly due to the very low 
luminosity of this source. 



5 CONCLUSIONS 

We have presented far-infrared Akari data of the L1155 re- 
gion at 90, 140 and 160 fim. We have compared these data 
to data of similar wavelengths as observed by the ISO satel- 
lite, and also to sub-mm data taken with SCUBA on the 
JCMT. The Akari data are shown to be consistent with the 
ISO data, although the Akari data cover a larger area. 

The data span the wavelength range of 90 - 850 fim, 
and have a monotonic shift in peak position of L1155C as a 
function of wavelength. We interpret this as a temperature 



gradient across the core. This hypothesis is backed up with 
SED fitting at the centre of the core and also at the front 
of the core, where the short-wavelength emission peaks. The 
best-fit temperature for the centre of L1155C is 11.9 K, com- 
pared to a temperature of 13.6 K at the front edge of the 
core. 

We carried out radiative transfer modelling of the 
L1155C core using the PHAETHON code. The code sub- 
jects a model core to the interstellar radiation field, plus a 
component from the nearby A6V star BD+67 1263. The core 
is treated as a plummer-like sphere, which has been skewed 
slightly towards the radiation source. With this model, we 
reproduce the appearance of the core, as observed at 90 - 
850 fim. We also reproduce both the temperature gradient, 
and the absolute temperature values. 

The LI 148 core exhibits a similar morphology at far- 
infrared and sub-mm wavelengths. Though we have fewer 
data available for this core, we are able to fit these data with 
a similar plummer-like sphere, which is externally heated 
from the nearby B-type star HD197053. 

The data presented here demonstrate that far-infrared 
data can very useful for interpreting longer wavelength data. 
In particular it can be useful for determining the temper- 
ature of cores, information which cannot be derived from 
sub-mm data alone. However, we note that far-infrared data 
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Figure 15. Graph of the emission from the model of L1148 at 
each observed wavelength, as a function of distance from the heat- 
ing source, measured along a straight line connecting the core and 
the heating source. The flux densities at both the centre and the 
front positions are plotted. 

should be used with caution for studying pre-stellar cores 
when complementary sub-mm data are not available. This 
study has shown how the density structure is not imme- 
diately apparent from the far-infrared data alone, and the 
sub-mm data are required to obtain a complete picture of 
the core. This needs to be borne in mind when studying 
data from other far-infrared satellites such as Spitzer and 
Herschel. 



ACKNOWLEDGMENTS 

The work is based on observations with AKARI, a J AX A 
project with the participation of ESA. 

The James Clerk Maxwell Telescope is operated by 
The Joint Astronomy Centre on behalf of the Science and 
Technology Facilities Council of the United Kingdom, the 
Netherlands Organisation for Scientific Research, and the 
National Research Council of Canada. The observations pre- 
sented here were taken under the programs M02AN18 and 
M03BN10. 

The Infrared Space Observatory is an ESA project with 
instruments funded by ESA member states particularly the 
PI countries, France, Germany, Netherlands and the United 
Kingdom and with the participation of ISAS (Japan) and 
NASA (USA). 

The Digitized Sky Survey was produced at the Space 
Telescope Science Institute under U.S. Government grant 
NAG W-2166. The images of these surveys are based on 
photographic data obtained using the Oschin Schmidt Tele- 
scope on Palomar Mountain and the UK Schmidt Telescope 
at Siding Spring. The plates were processed into the present 
compressed digital form at the Royal Edinburgh Observa- 
tory (ROE) photolabs with the permission of these institu- 
tions. 

The authors acknowledge the use of NASA's SkyView 
facility (http://skyview.gsfc.nasa.gov) located at NASA 
Goddard Space Flight Center. 



This research made use of Montage, funded by the Na- 
tional Aeronautics and Space Administration's Earth Sci- 
ence Technology Office, Computation Technologies Project, 
under Cooperative Agreement Number NCC5-626 between 
NASA and the California Institute of Technology. Montage 
is maintained by the NASA/IPAC Infrared Science Archive. 

The authors wish to thank J. Kirk for the use of his 
SED fitting scripts, and Nicholas Chapman for the use of 
his Spitzer-IRAC extinction map of L1155C. 

DN and DS acknowledge STFC for PDRA support un- 
der the Cardiff University Astronomy Rolling Grant. 



REFERENCES 

Andre P., Ward-Thompson D., Motte F., 1996, A&A, 314, 
625 

Andre P., Ward- Thompson D., Barsony M., 2000, 
prpl.conf, 59 

Andre, P.,Bouwman, J., Belloche, A., & Hennebelle, P. 

2003, in Chemistry as a Diagnostic of Star Formation, 
University of Waterloo, Canada. Eds. Charles L. Curry 
and Michel Fich. NRC Press, p. 127. 

Arce H. G, Goodman A. A., 2002, ApJ, 575, 911 

Archibald E. N., et al, 2002, MNRAS, 336,1 

Bacmann A., Andre P., Puget J.-L., Abergel A., Bontemps 

S., Ward-Thompson D., 2000, A&A, 361, 555 
Baes M., Stamatellos D, Davies J. I., Whitworth A. P., 

Sabatini S., Roberts S., Linder S. M., Evans R., 2005, 

New A, 10, 523 

Balazs L. G, Abraham P., Kun M., Kelemen J., Toth L. V., 

2004, A&A, 425, 133 

Bjorkman J. E., Wood K., 2001, ApJ, 554, 615 
Black J. H., 1994, ASPC, 58, 355 

Dobashi K., Uehara H., Kandori R., Sakurai T., Kaiden 

M., Umemoto T., Sato F., 2005, PASJ, 57, 1 
Emerson D. T., 1995, in ASP Conf. Ser. 75, Multifeed Sys- 
tems for Radio Telescopes, ed. D. T. Emerson & J.M. 
Payne (San Francisco: ASP), 309 
Goodman A. A., Arce H. G, 2004, ApJ, 608, 831 
Griffin M. J., Orton G. S., 1993, Icar, 105, 537 
Hildebrand R. H., 1983, QJRAS, 24, 267 
Holland W. S., et al., 1999, MNRAS, 303, 659 
Jenness T., Lightfoot J.F., 2000, Starlink User Note 216, 

Star link Project, CCLRC 
Jessop N. E., Ward- Thompson D, 2001, MNRAS, 323, 
1025 

Kawada M., et al., 2007, PASJ, 59, 389 
Kessler M. F., et al., 1996, A&A, 315, L27 
Kirk J. M., Ward-Thompson D., Andre P., 2005, MNRAS, 
360, 1506 

Kirk J. M., Ward-Thompson D., Andre P., 2007, MNRAS, 

375, 843 
Kirk J. M. et al., 2008, in prep. 

Klaas U., Kruger H., Heinrichsen I., Heske A., Laureijs 
R. J., 1994, ISOPHOT Observers Manual. ESA, Noord- 
wijk 

Kun M., 1998, ApJS, 115, 59 

Kun M., Prusti T., 1993, A&A, 272, 235 

Kun M., Vinko J., Szabados L., 2000, MNRAS, 319, 777 

Lasker B. M., 1994, IAUS, 161, 167 



Akari mapping of pre- stellar cores 15 

Li W., Evans N. J., II, Harvey P. M., Colome C, 1994, 
ApJ, 433, 199 

Looney L. W., Tobin J. J., Kwon W., 2007, ApJ, 670, L131 
Lynds B. T., 1962, ApJS, 7, 1 

Kauffmann J., Bertoldi F., Evans N. J., II, the C2D Col- 
laboration, 2005, AN, 326, 878 

Kun M., Kiss Z. T., Balog Z., 2008, arXiv, arXiv:0809.4761 

McGlynn T., Scollick K., 1994, ASPC, 61, 34 

Murakami H., et al., 2007, PASJ, 59, 369 

Myers P. C, Linke R. A., Benson P. J., 1983, ApJ, 264, 
517 

Neugebauer G., et al., 1984, ApJ, 278, LI 
Nutter D., Ward-Thompson D., 2007, MNRAS, 374, 1413 
Motte F., Andre P., Neri R., 1998, A&A, 336, 150 
Orton G. S., Griffin M. J., Ade P. A. R., Nolt I. G., Ra- 

dostitz J. V., 1986, Icar, 67, 289 
Perryman M. A. C, et al., 1997, A&A, 323, L49 
Plummer H. C, 1911, MNRAS, 71, 460 
Racine R., 1968, AJ, 73, 233 
Sandell G., 1994, MNRAS, 271, 75 

Simpson R. J., Nutter D., Ward- Thompson D., 2008, MN- 
RAS, 391, 205 
Snell R. L., 1981, ApJS, 45, 121 

Stamatellos D., Whitworth A. P., 2003, A&A, 407, 941 
Stamatellos D., Whitworth A. P., Ward-Thompson D., 

2007, MNRAS, 379, 1390 
Straizys V., Cernis K., Kazlauskas A., Meistas E., 1992, 

BaltA, 1, 149 
Tafalla M., Bachiller R., 1995, ApJ, 443, L37 
van den Ancker M. E., The P. S., Tjin A Djie H. R. E., 

Catala C, de Winter D., Blondel P. F. G, Waters 

L. B. F. M., 1997, A&A, 324, L33 
Verdugo E., Yamamura I., Pearson C, 2007, Akari FIS 

Data User Manual vl.3, http://www.ir.isas.ac.jp/ 
Viotti R., 1969, MmSAI, 40, 75 

Ward- Thompson D., Scott P. F., Hills R. E., Andre P., 

1994, MNRAS, 268, 276 
Ward- Thompson D., Motte F., Andre P., 1999, MNRAS, 

305, 143 

Ward- Thompson D., Andre P., Kirk J. M., 2002, MNRAS, 
329, 257 

Ward- Thompson D., Andre P., Crutcher R., Johnstone D., 

Onishi T., Wilson C, 2007, prpl.conf, 33 
Whitworth A. P., Ward- Thompson D., 2001, ApJ, 547, 317 
Yonekura Y., Dobashi K., Mizuno A., Ogawa H., Fukui Y., 

1997, ApJS, 110, 21 



